Refinement of the equation of state of tantalum 
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The volume of tantalum versus pressure has been accurately measured up to 101 GPa by single- 
crystal x-ray diffraction, with helium as pressure transmitting medium. Slight deviation from pre- 
vious static determinations is observed. Discrepancy with reduced shock-wave and ultrasonic data 
supports recent doubts about the calibration of the ruby pressure scale. Finally, first principle cal- 
culations of the literature show a positive curvature in P(V) relative to the experimental data, even 
with a modified pressure scale. 
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Ta is a transition metal (Z=73 and the Xe 4f 14 5d 3 
6s 2 equilibrium atomic configuration), with a simple bcc 
structure stable to at least 174 GPai. Because of its 
fundamental interest and its use as an important tech- 
nology material, the equation of state (EoS) of Ta has 
been studied by several experimental techniques 2 ^ 3 ")^^ 
and numerous ab-initio density functional theory (DFT) 
electronic structure calculations^!^. However, the EoS 
of Ta is still not known with a satisfactory accuracy. On 
the experimental side, small differences between static 
determinations are observed^, and more important dis- 
crepancies exist among ultrasonic measurements 2 -, shock- 
wave results 3 -^, and static determinations^. On the 
theoretical side, systematic deviations are seen between 
calculations with the two mostly used approximations of 
the exchange-correlation energy of the electrons, and re- 
liable experimental data are needed to validate one of 
these approximations^. Here below, we report accurate 
EoS data of Ta, based on single-crystal x-ray diffraction 
(XRD), up to 101 GPa under quasi-hydrostatic condi- 
tions. We compare these data with previous measure- 
ments and DFT calculations. 

Six runs were dedicated to the measurement of accu- 
rate volume data of Ta versus pressure by synchrotron 
single crystal x-ray diffraction, on the ID30 beamline at 
the ESRF (Grenoble, France). Small single crystal grains 
(4 /im in the maximum dimension) , chosen in a tantalum 
powder (averaged size 10 /im, 99.9 % purity, Goodfel- 
low product) on the basis of their external shape, were 
loaded in membrane diamond anvil cells, with helium as 
pressure transmitting medium. We checked by interfer- 
ometry that the thickness of sample chamber was always 
larger than the dimension of the crystal. The first two 
experiments have been carried out using single crystal 
energy dispersive XRD and the other four using angle 
dispersive monochromatic XRD at 0.3738 A. The same 
technique had been used to measure with high accuracy 
the equation of state of low-Z systems in the 100 GPa 
rangeiS. An average of 7 reflections were measured in 
each run ; absolute uncertainty in the lattice parame- 
ter is at maximum of 10 -3 . The pressure was estimated 
from the luminescence of a small ruby ball (less than 4 
fim in diameter) and its quasi-hydrostatic calibration^. 



Non-hydrostatic stresses were observed to be negligible 
because the shape of the diffraction spot for a given re- 
flection showed no sign of deterioration of the crystalline 
quality of the sample. Except in one run, when a leak 
of helium through a diamond fracture caused a deteri- 
oration of the pressurization conditions above 70 GPa, 
and the volume of a strained sample, bridged between 
the diamond tips, was subsequently measured. 
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TABLE I: Atomic volume of Ta measured versus pressure by 
single-crystal XRD in helium pressure transmitting medium. 



The P — V data points measured under hydrostatic 
conditions are presented in Table U They are plotted in 
Fig. and compared to the data obtained with uniax- 
ial stress conditions in the broken diamond run. It is 
seen that a non-hydrostatic pressure leads to overesti- 
mate the volume. The ambient pressure atomic volume 
was measured equal to 18.033(50) A 3 , which agrees with 
ambient T and P literature data (18.065 A 3 ,12). The 
uncertainty in the volume determination is smaller than 
3.10 -3 , which corresponds to a pressure uncertainty of 
1.5 GPa at 100 GPa. The estimated pressure error bars 
include the pressure gradient between the ruby ball and 
the Ta crystal and range from 0.05 GPa at 1 GPa to 1 
GPa at 100 GPa, if the ruby pressure scale is assumed 
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FIG. 1: Pressure versus volume data of Ta at 300K. The 
dots and circles are respectively the data measured under hy- 
drostatic conditions and non-hydrostatic conditions. The full 
line is the Vinet fit of the data as discussed in the text. P = 
3K x- 2 (l-x)exp(1.5(K[ l -l)(l-x)), with x = 1-(V/Vb) 1/3 
and Vo = 18.033 A 3 , Ko = 194 GPa, K' Q = 3.25. Inset : 
F(x) = \n(Px 2 /(l - x)) = ]n(3K ) + %(K - - x), is 
the linearized form of the Vinet function which is used to 
determine K' ±0.1. 
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FIG. 2: Pressure difference between various experimental de- 
termination of the pressure-volume relationship of Ta. AP 
expresses the difference for a given atomic volume, between 
the various experimental data and our fitted Vinet EoS, and 
it is plotted versus the atomic volume (bottom) and the cor- 
responding pressure given by the Vinet EoS discussed in the 
text. 



tion of the ruby pressure scaled) are compared in ta.hlelTTl 
to ultrasonic measurements (the adiabatic to isothermal 
correction had been done in refi^) and to the reduced 
Hugoniot datal. The Kq value obtained in the current 
study, 198 ± 3.4 GPa, is in correct agreement with the 
ultrasonic one, 194 GPa. On the contrary, a smaller Kq 
is obtained, wether Ko is fixed or not to its ultrasonic 
value. 



to be correct. Therefore, the error in the volume de- 
termination, intrinsical to the XRD technique is mainly 
responsible of the scattering of the experimental points 
as seen in Fig. 

To provide the useful physical parameters, namely vol- 
ume Vo, bulk modulus Kq, and its pressure derivative 
K Q , under ambient conditions, the P — V data points are 
generally fitted by an EoS formulation. The least square 
fit of the data with an universal form of the equation of 
state should ideally give the same value for Vo, Ko and 
Kq as the ones measured independently at zero pressure 
by XRD for Vo and by ultrasound propagation for Kq 
and K . We have fitted the data with three most-used 
forms of the EoS in the literature: third-order Birch- 
Murnaghan EoSH (BM), Vinet EoSl4 (V) and the EoS 
proposed by Holzapfeli^ (H). In this pressure range, very 
similar values of K (between 196 and 199 GPa), and K Q 
(between 3.07 and 3.30), are obtained with these three 
forms, that give comparable quality of the fits. The Vinet 
fit has been chosen to reproduce the present data. In the 
case of the Vinet fit, the data can be seen in a linearized 
form, as shown in the inset of Fig. It is seen that the 
data points above 50 GPa are important to unambigu- 
ously constrain the Kq value of the slope. The results of 
the various Vinet fits of the data (fixing Vo, as measured; 
fixing Vo and K as given by ultrasonic measurements; 
using the recently proposed modification of the calibra- 
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TABLE II: Values of V (in A 3 ), K (in GPa) and K' Q for 
the EoS of Ta. The current and literature X-Ray diffraction 
(XRD) measurements are compared to ultrasonic (US) and 
reduced shock-wave (RSW) measurements. For each mea- 
surement, the technique, the pressure range and the form of 
the EoS (V, BM, H) are indicated. The parameters that have 
been fixed in the fitting procedure are in bold font. Numbers 
between parenthesis are the fitting error bars (95 % confidence 
interval). The last line presents the EoS parameters obtained 
after modification of the pressure scaled. 

In Fig. other experimental determinations of the 
EoS of Ta are compared through their differences with 
the Vinet fit of the present data (with Vb=18.033 A 3 , 
Kq = 194 GPa, Kq = 3.25). A systematic deviation be- 
tween the present data and the ones of two recent XRD 
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measurements is clearly seen. The numerous low pres- 
sure measurements of Hanfland et al. give an atomic vol- 
ume at ambient pressure of 18.186 A 3 (0.8 % larger than 
the present determination) . This could be due to a small 
contamination by hydrogen (a few percents) , that is well- 
known to lead to an expansion of the lattice^. For this 
reason, the points plotted in Fig. [51 appear shifted from 
our data. Due to its larger compressibility at low pres- 
sure (Fig. |2J , the values of Kq fitted on Hanfland data 
cannot be reconciled with the ultrasonic value. However, 
the compression curve of Hanfland et al. is parallel to 
the current one at higher pressure. Low pressure data 
points measured by Cynn and Yooi exhibit a large scat- 
ter. Also, an increasing deviation between the present 
EoS data and the data of Cynn is observed at P > 60 
GPa. This could be ascribed to an uniaxial stress effect 
in Cynn experiments^ because these data are following 
the trend of our data points in the non-hydrostatic run 
(in this study, gold has been used as an x-ray pressure 
standard^ for the high pressure data points) . 

More severe inconsistencies exist between the present 
XRD and the ultrasonic measurements or the reduced 
shock-wave EoS of tantalum. All XRD studies lead to 
a smaller Kq than the ultrasonic one. Inserting the ul- 
trasonic values^ of Kq and Kq in the Vinet formulation 
generates an equation of state that deviates from our de- 
termination at high pressure. As seen in Fig. 2, the EoS 
generated from ultrasonic values of Kq and K is in very 
good agreement with the EoS of Ta obtained by the re- 
duction of shock data to ambient UljA. The difference 
between the XRD EoS and the reduced shock EoS or 
the ultrasonic based EoS reaches 12 % in pressure at 150 
GPa. Possible causes of these apparent inconsistencies 
need to be examined and error bars estimated. The er- 
ror on the reduced Hugoniot data of tantalum^ is claimed 
to be of 5.5 % at 400 GPa, and less below 200 GPa. This 
uncertainty is caused both by uncertainties in particle 
and shock velocities measurements, and by the reduction 
procedure from Hugoniot temperature to the 300 K. The 
sum of the error bars in the volume determination and in 
the pressure measurements on the XRD EoS data points 
amounts to an equivalent 2.5 GPa at maximum at 100 
GPa. Consequently, the consideration of cumulative er- 
ror bars on the reduced Hugoniot and on the XRD EoS 
cannot reconcile the dynamic and static determinations 
of the EoS of Ta. To reconcile our XRD data with the 
reduced Hugoniot data and the ultrasonic based EoS of 
tantalum, it must be assumed that the ruby scale un- 
derestimates the pressure by about 10% at 150 GPa. A 
similar correction of the ruby pressure scale has been re- 
cently proposed by Holzapfe!~. Moreover, a discrepancy 
between ab-initio calculations and x-ray measurements 
of the pressure- volume relation of diamond has also been 
ascribed to an error of the ruby pressure scale of the 
same magnitude^. As seen in table HH by using the new 
ruby pressure scale proposed by Holzapfel, consistency 
is recovered between the EoS of Ta obtained by static 
XRD and shock-wave measurements because similar val- 
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FIG. 3: Comparison between DFT calculations and experi- 
mental P-V data points at 300K. The pressure difference, for 
a given atomic volume, between the ab-initio pressure and 
the pressure of our fitted Vinet EoS is plotted versus volume. 
GGA and LDA calculations^ are respectively represented by 
bold and thin lines. The pressure change due to ruby scale 
scale of Holzapfel— is indicated as dotted line 



ues of Kq, and K' are then obtained. The ultrasonic 
K value remains slightly (0.28) larger than the one ob- 
tained by fitting of compression data, which could be due 
to either experimental uncertainties or approximations in 
finite strain EoS formulation. 

The ab-initio calculations of the EoS of Ta has been 
produced in four recent studiesi^^^ using density- 
functional theory. In principle, this method involves ap- 
proximations only in the exchange and correlation en- 
ergy calculation. Two forms are currently used, namely 
the Local Density Approximation (LDA) and the Gener- 
alized Gradient Approximation (GGA). The validity of 
these approximations is in practice established by their 
ability to reproduce experimental results. Uncertainties 
might also arise consequently to the choice of the com- 
putational method and of its inherent assumptions, but 
they should be small, as it is seen by comparing three 
recent GGA EoSi^. These EoS agree within 3 GPa 
at 150 GPa, although they are based on different com- 
putational methods. Consequently, for legibility in Fig. 
3, the LDA and GGA EoS are represented by the cal- 
culation of Boettgei4i alone, and they are compared to 
the experimental data through the pressure difference at 
a given volume. It is observed here, as often, that the 
ambient pressure volume is underestimated by LDA cal- 
culations and overestimated by GGA calculations. What 
seems to be more interesting here is that all DFT calcu- 
lations show a similar positive curvature relative to the 
experimental data. The higher compressibility at high 
pressure of the experimental EoS compared to the DFT 
calculations cannot be entirely due to the possible error 
of the ruby pressure scale. In fact, if the new calibration 
proposed by Holzapfel^, that reconciles all the experi- 
mental determinations as discussed above, is used, the 
larger experimental compressibility still exists, as shown 
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in Fig. for pressures higher than 80 GPa. 

In summary, the present work demonstrates that an 
accurate determination of the volume of metals under 
hydrostatic conditions in the 100 GPa range can now be 
achieved by synchrotron single-crystal XRD with helium 
as pressure transmitting medium. Consistency with the 
determination of the EoS of Ta by ultrasonic or shock- 
wave methods can be recovered by using a calibration 
of the ruby scale recently proposed. This highlights and 
estimates the present uncertainty of the ruby pressure 
calibration. Similar EoS measurements under good hy- 
drostatic conditions up to the Mbar range should now 
be performed on reference metals to address two impor- 
tant issues pointed out in the present work: first, to re- 
duce the uncertainty of the ruby pressure scale, by a 
comparison of these metal EoS data with the EoS de- 
duced from high-temperature shock-wave, considered as 



primary standardsii; second, to investigate if a higher ex- 
perimental compressibility than the calculated DFT one 
is a systematic trend in metals and if it has a Z depen- 
dence. 
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